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Within a spin wave analysis and a fermionic large-n limit, it is shown that the antiferromagnetic 
Heisenberg model on the checkerboard lattice may have different ground states, depending on the 
spin size S. Through an additional exchange interaction that corresponds to an inter-tetrahedra 
coupling, the stability of the Neel state has been explored for all cases from the square lattice to 
the isotropic checkerboard lattice. Away from the isotropic limit and within the linear spin wave 
approximation, it is shown that there exists a critical coupling for which the local magnetization of 
the Neel state vanishes for any value of the spin S. One the other hand, using the Dyson-Maleev 
approximation, this result is valid only in the case S = ^ and the limit between a stable and an 
unstable Neel state is at 5 = f. For S — ^, the fermionic large-n limit suggests that the ground 
state is a valence bond solid build with disconnected 4-spins singlets. This analysis indicates that for 
low spin and in the isotropic limit, the checkerboard antiferromagnet may be close to an instability 
between an ordered 5 = ground state and a magnetized ground state. 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.50.Ee 



I. INTRODUCTION 

In one dimension (ID) the ground state of the Heisen- 
berg antiferromagnet may be quasi-ordered or disor- 
dered, depending on the spin parity SEJ. When S is an 
integer, and in particular, when S* = 1, the ground state 
corresponds to a magnetically disordered state where 
spin-spin correlation functions (CFs) exponentially de- 
cay with distance and with total spin St — 0- This state 
is called a quantum spin liquid (QSL) as the spin-spin 
CFs are equivalent to the density-density CFs in conven- 
tional liquids and because it has total spin S = and 
no broken symmetry. This concept can be generalized to 
classical systems if the spins are considered as 0(3) vec- 
tors. The constraint 5 = is then removed and the state 
is termed a classical spin liquid (CSL). Another notion for 
spin liquid ground states emerged from earlier studies by 
Wanniercl and Houtappelu and is related to geometrical 
frustration. In that case, it is the geometry of the un- 
derlying lattice that leads to magnetic disorder and not 
its dimensionality. Naturally, frustrated structures may 
exist not only for D = 1 and for over 50 years frustrated 
lattices have been examined in the search for disordered 
ground states in higher dimensions {D = 2, 3). One goal 
of this research into peculiar structures is the study of 
unconventionnal spectra, not covered by the Lieb and 
Mattis theorenJj. 

To date, two particular lattices have been considered 
as good candidates for spin liquid ground states : the 
kagome lattice and the pyrochlore lattice. The first one 
is a two dimensional arrangement of corner sharing trian- 
gles, while the latter is a three dimensional structure of 
corner sharing tetrahedra. Both lattices have been the- 
oretically investigated and it has been suggested that in 
the classical well as in the quaBtiim case, they 

possess spin liquid like ground statesBEiCHa. Experimen- 



tally, all systems that can be well described by the an- 
tiferromagnetic Heisenberg model on these lattices have 
been found to possess this spin liquid behavior or un- 
convential behaviors that can be ascribed in part, to the 
underlying "theoretical" spin liquid ground stateQ. This 
suggests that both the elementary cells (triangle or tetra- 
hedron) and the connectivity (corner sharing) are ma- 
jor ingredients for magnetic disorder. In this paper, we 
study such a lattice, where the elementary cell is a tetra- 
hedron and where connectivity is corner sharing : the 
checkerboard lattice (see Fig. p. This system is two di- 
mensionnal and can be described as the two dimensional 
analog of the pyrochlore lattice : each atom possess lo- 
cally the same environment. It has recently been proven 
that its quantum ground states are singletslij for finite 
systems. The aim of this work is to test wether this con- 
clusion may be valid for an infinite lattice. To do so we 
start away from the isotropic checkerboard limit, when 
diagonal bonds are weakened (see Fig. and we test 
wether this system undergoes a Neel ordering transition 
when reaching the isotropic limit. This approach is very 
close irt^irit to the one previously done by Chandra and 
DouQoilB on the J1-J2 square lattice model. This is done 
within the frame of the linear-spin wave theory and the 
Dyson-Maleev approximationllj. Whereas for the linear 
spin wave approach there always exists a critical coupling 
for which the local magnetization is unstable, whatever 
the value of the spin, in the case of the Dyson-Maleev 
approach there is a critical coupling only for the case 
S — ^. The peculiar case S" = ^ is explored within a 
fermionic largc-n limit of the SU(2) hamiltonian. It is 
found that the system should order in a valence bond 
solid of disconnected 4-spins singlets. The notations, the 
hamiltonian, and the linear spin wave analysis are intro- 
duced in Sec. ||. The Dyson-Maleev formalism is shortly 
described in Sec. as it is can be seen as an extension 
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FIG. 1: The checkerboard lattice. It can be pictured as a 
square lattice of tetraheda (top). Locally, it reproduces the 
same environment as the three dimensional pyrochlore lattice. 
In order to simplify our calculations, it has been described by 
a square lattice with a 2-spins unit cell (bottom). 



of Sec. |l^ Results of the linear spin wave theory and the 
Dyson-Maleev approach are presented and discussed in 
Sec. [V, The case 5 = ^ is detailed in Sec. ^. 



II. LATTICE, HAMILTONIAN AND SPIN 
WAVE ANALYSIS 



A. Notations and Lattice Description 

The checkerboard lattice can be described as a square 
lattice with a unit cell containing 2 spins (see Fig. |^). 
Each cell of the lattice can be obtained from the others 
by applying the translations a = (1,0) and b = (0,1). 
This means that each site of the lattice can be defined by 
two indices {i,n). The first corresponds to the unit cell 
(i = 1 . . . N) where N is the total number of cells in the 
lattice and the other to the type of the site {n — 1,2). 
We introduce two kinds of interactions, J and J'. J cor- 
responds to a coupling constant on a square lattice. J' 
can be seen as a coupling constant within one square over 
two of the square lattice (see Fig. ||). Using these nota- 
tions, the Brillouin zone of the underlying Bravais lattice 
is BZ(kx, ky) = [— TT, tt] X [— tt, tt] and the structure coeffi- 
cients can be defined a.s 5 = cos{k.S) — cos ((fcx + fcy)/2). 



FIG. 2: (left) Lattice of coupling constants J and J'. Bold 
lines design constant J while dotted lines describe constant 
J', (right) Local environment of a unit cell. 5 and 7 are 
related to J- like coupling while e* and ff correspond to J'. 



7 = cos(/c.7) 
rj = cos(fc.rf) 
in Fig. y. 



y)/2), 



COS ((/cx + k 
cos(fcy) where 6 



— cos{k.e) — cos{k^), 
e and ff are depicted 



B. Spin Wave Hamiltonian 

The Heisenberg hamiltonian for this system is 

H ~ —"^^JijSi.Sj (1) 

ij 

where Jij are coupling constants {Jij < for an antifer- 
romagnet) and Si is the spin located at site i. This spin 
hamiltonian may be transformed into a boson hamilto- 
nian through an Holstein-Primakoff transformation, pro- 
vided that we know its Mean Field ground states (MF) 
from the square lattice limit to the checkerboard lattice 
limit, i.e for all values of r = J'/ J. One way to know the 
MF ground state(s) is to compute, the Fourier transform 
of the interactions on the latticet^l. 



H 



ijSi.Sj 



^ J(g)%S'_,- 



(2) 



As we here deal with a 2-spin unit cell, its Fourier trans- 
form J{q) is a 2 X 2 matrix and must be diagonalized to 
find the normal modes X^^^'^{q). The ground state (if 
unique) then corresponds to the lowest eigenvalue X^{qo) 
and is described by the propagation vector qo and the 
associated eigenvector Vi{qo) which fixes the inner struc- 
ture of the unit cell. In the range < r < 1, (fo always 
equals and the eigenvector is always equal to (1,-1) 
which means that the ground state is non-degenerate 
and corresponds to the antiferromagnetic Neel state (see 
Fig. H top and middle). At the peculiar point r = 1, we 
obtain a flat branch, signaling that at the MF level, the 
system is continuously degenerate (see Fig. || bottom), 
as previously noted in Ref. |l^. That behavior can be 
explained in terms of local degrees of freedom. In the 
limit J' — J, we note that the hamiltonian may be de- 
composed as a sum of squares over "tetrahedra" , as it 
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FIG. 3: Fourier transform of the interactions for the square 
lattice limit (r — 0, top), an intermediate case (r = 0.5, 
middle) and the checkerboard limit (r = 1, bottom). Each 
manifold corresponds to an eigenvalue of J{q). 

has been done in the kagome as well as in the pyrochlore 
lattice 

H = — J ^ " ^ Sj-Sj 

ij 

= -?EfE^'l +2NJS(S + 1) 



(3) 



Tot. \i=l 



FIG. 4: Fourier transform of the interactions for r = 2. Each 
manifold corresponds to an eigenvalue of J{q)- The mean 
field ground state is continuously degenerate along lines which 
unable the choice of a trial ground state for the spin wave 
analysis. 



We see that the classical ground state is defined by the 
constraint S'i = on each tetrahedron. Further- 

more, as tetrahedra are corner sharing, this constraint is 
not sufficient to order the system and there are an infi- 
nite number of ways to satisfy it, which results in a flat 
branch over the whole Brillouin zone. This is exactly the 
same phenomenon that was previously observed on the 
kagome as well as on the pyrochlore lattice. Nevertheless, 
the Neel state belongs to this infinite set of MF ground 
states, and as it is the correct one for any r < 1, we 
will also consider this state by continuity for the point 
r = 1. For the parameter range r > 1, the situation is 
quite different as MF ground states are continuously de- 
generate along lines corresponding to the border of the 
square lattice Brillouin zone (see Fig. ^). This precludes 
a description via a spin wave analysis in this range as no 
magnetic ground state can be choosen as a trial state. 
This is the same situation as the one encountered on the 
kagome or the pyrochlore lattices where spin wave calcu- 
lations can only compare relative stability of NeeL-states 
without considering the entire set of these statesEJ. So, 
we will restrict ourselves to the range < r < 1 , and con- 
sider the antifcrromagnetic Neel state as the trial state. 
We therefore introduce two species of bosons, each one 
related to a different sublattice, and restrict the devel- 
oppement of spin operators to the lowest order in the 
spin wave expansion The resulting bosonic hamiltonian 
is then diagonalizcd through a Bogoliubov transforma- 
tion. We are left with two branches of "magnons" whose 
dispertions are = (l/2)(/ ~ g + y/{f -I- gY - AK^) 
and ep = (l/2)(-(/ - g) + ^{f + gY-Ah^) where 
/fe = 2 + r(e - 1), .gfc = 2 + r(r? - I), hk = 5 + i and 
r = J' I J ((5, 7, e and r] are defined in Sec. II A ). 

Restricting to the ground state, the energy as well as 



the local magnetization are obtained through 

E = = 2JS^ - J'S^ - JS'^Cik) (4) 

k 

M = {St)^^ S - {4c,)^ = S - ^Y.'^lik^ (5) 

fc 

with Ck = i(ea + e/3 + (/ - g)iul ~ uj) - (/ + 5)) 
and where the coefficient s of the Bogo liubov trans- 
formation are Ua = h/ ^h? — (eq — /)^, Va = if — 

and = h/^{ep + f)^-h^. 



III. DYSON-MALEEV REPRESENTATION 

We use the same hamiltonian of Kq. |l] and intro- 
duce the Dyson-Maleev transformationtj. Because of its 
algebra, the Dyson-Maleev transformation rewrites the 
hamiltonian of Eq. |^ as a sum of two and four opera- 
tors terms. The simplest method to solve the problem 
is to decouple the four operators term introducing an a 
priori decoupling scheme. In the present case, we follow 
the mean field decoupling of Ref. Three quantities 
are introduced that correspond each to a two bosons like 
excitations, 

V = (a+a,) = {bp^) 

O = (a+a,) = {b+bj ■ , j^m (6) 

Other types of excitations are neglected and consid- 
ered as higheiL-cnergy processes as previously proposed 
by TakahashiEj, Using these approximations, the Dyson 
Maleev hamiltonian is then linearized and diagonalized 
through another Bogoliubov transformation. The struc- 
ture of the transformation is identical to the one of the 
linear spin wave case, but with new coefficients. These 
coefficients are fl = 2A + rB{e-l), = 2A + rB{T]-l), 
h',^ = A{S + 7) where A = 1 - (V + 0)/2S and B = 
1 — {T> — T)/2S. Finally the three self consistent coeffi- 
cients at zero temperature are obtained through the three 
coupled relations, 

k 

k 

k 

and the local magnetization of the Neel state is simply 
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FIG. 5: Energy per site as a function of the ratio r = J' / J, i.e 
between the limits of the square lattice and the checkerboard 
lattice, for S = 1/2. Beyond the critical ratio r « 0.75, the 
value of the energy is meaningless as it describes the energy 
of an unstable state. 
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FIG. 6: Local magnetization as a function of the ratio r = 
J'/ J, i.e from the square lattice to the checkerboard lattice. 
Red lines correspond to the calculations made for 5=1/2 
and black lines for 5 = 1. 



IV. RESULTS OF THE SPIN WAVE 
APPROACH 

In the following, we fix J = —1 and take — 1 < J' < 0. 
We first discuss the linear spin wave results and then the 
Dyson-Maleev approach. In Fig. || and ^, the energy as 
well as the local magnetization are plotted as functions 
of the ratio r — J' /J. For the case r = 0, well known 
results for the S ~ 1/2 square lattice are recovered : 
E — —0.658 and AI — 0.3. The magnon spectrum is very 
simple; there is only one dispersive branch, which is two 
times degenerate in our case because of the 2-spin unit 
cell description (Fig. ^top). With increasing r, there are 
three zones to distinguish. First, in the range < r < Tc, 
the magnon spectrum has two distinct manifolds, each of 
these being dispersive (Fig. ^ middle). This occurs now 
because the site of the unit cell no longer have the same 




FIG. 7: Magnons dispersions ecjiik) over the Brillouin zone. 
From top to bottom, r = 0, 0.5, 1. At r = 1, the lowest branch 
is non-dispersive. This means that this point is remarkable 
and reveals (as in the kagome and the pyrochlore lattice case) 
the presence of a soft mode related to local degrees of freedom. 
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FIG. 8: Stability of the Neel state. The straight line indicates 
the threshold below wich the Neel state is stable for a given 
ratio r — J' I J . Over this line, it is destabilized by geometry 
and quantum fluctuations. For r = 1, it is always unstable, 
whatever the value of the spin S. 



symmetry. Tc corresponds to a critical coupling below 
which the local magnetization is finite, even renormahzed 
by quantum corrections. For S = 1/2, « 0.75 which 
was already obtained in Ref. |l^ and for S = 1, Tc ~ 0.93. 
Critical values for higher spins are shown in the phase 
diagram of Fig. |^. The second zone corresponds to the 
range Tc < r < 1. Because of quantum corrections, local 
magnetization is renormalized to zero, which means that 
in this approximation, the system does not show mag- 
netic ordering. Even if quantum fluctuations drive the 
system to disorder, one notes that magnons are still dis- 
persive. This means that there is no peculiar contribution 
coming from their structure. This is no longer the case for 
the special point r = 1. At this point, the lowest branch 
of spin waves is flat over the whole Brillouin zone (Fig. 
bottom) . This has several consequences. First, this 
means that at the level of this approximation, there will 
always be a critical coupling, whatever the value of the 
spin. This comes directly from the non-analytic behavior 
olva{k) when r ^ 1 along directions ky = ±kx. Because 
of that the local magnetization diverges as — r), 

i.e is renormalized to zero, whatever the value of S. Con- 
sequently, we obtain a phase diagram {1/S,J' / J) very 
similar to the one of the Ji — J2 square lattice model 
where J'/ J plays the role of J2/J1 (see Fig. || and Fig. 2 
of Ref. O). This problem of flat mode is often encoun- 
tered in frustrated lattices and it rules out the possibility 
of studying the r = 1 point without going to next order 
of expansion in our model, i.e, without including quar- 
tic terms. Whether such an inclusion of quartic terms is 
relevant is not clear as it corresponds to consider quartic 
modes as perturbations of zero quadratic modes. Nev- 
ertheless, it is the simplest next step to be done-in the 
frame of spin wave (Holstein-Primakoff) analysisEJ. Sec- 
ond, it has a meaning very close to MF conclusions. If 
the lowest branch of magnons is flat, this means that the 



FIG. 9: Relative local magnetization of the Neel state ver- 
sus the inverse spin value 1/5 at the isotropic point. For all 
S > 1, the local magnetization is finite, suggesting that the 
ground state could be the Neel state. For S — 1/2, the local 
magnetization is renormalized to zero by quantum fluctua- 
tions. 



system has no spin stiffness. Thus, any local deformation 
will not propagate through lattice. This is very close to 
the argument of continuous degeneracy obtained through 
the MF approach, i.e the existence of an infinite number 
of local degrees of freedom (see Eq. ||). These degrees 
of freedom are one of the signatures of frustrated sys- 
tems and may be counted, in the case of classical spins, 
through the "Maxwellian counting" ruleQ. At the level of 
the LSW approximation, this means that quantum fluc- 
tuations and geometrical frustration may destabilize the 
Neel state for a finite range of J'. Wether a quantum 
spin liquid ground state in this range is stable is beyond 
this approach even if the existence of critical couplings 
goes in the right direction. 

We now turn to the Dyson-Maleev approach. As for 
the preceding case, the local magnetization has been com- 
puted for different value of the spin. The main result is 
that for all S' > 1, there is no critical coupling, i.e the 
local magnetization of the Neel state is always finite in 
the isotropic limit r = 1. This is illustratred in Fig. ^ 
where reduced local magnetization {Sf)f^/S is plotted 
versus the inverse spin value 1/S at r —>■ 1 . The only 
exception is for S ~ ^ for which it seems that there is 
a critical coupling at r « 0.98. To test wether the sta- 
bility of the Neel state for the S = ^ comes from a fiat 
mode or from quantum fluctuations, we have plotted the 
associated dispersions over the whole Brillouin zone (see 
Fig. |l^). From the structure of the bosons excitations, 
it is clear that including quantum fluctuations via the 
Dyson-Maleev formalism destroys the pathological struc- 
ture of the flat mode. Therefore, destabilization of the 
Neel state is only driven by quantum fluctuations within 
this approximation and not by a zero spin stiffness as 



FIG. 10: Energy dispersions corresponding to the bosons ex- 
citations around the Neel state for the S — 1/2 case in the 
limit r ^ 1. The dispersions are no longer pathological as 
they are in the LSW approximation and the flat mode is de- 
stroyed. This clearly indicates that destabilization of the Neel 
state is driven here only by quantum fluctuations. 



it was the case in the LSW approximation in the limit 

These results are quite surprising as they are quali- 
tatively different from the LSW approximation. This 
suggests that the more accurately the quantum fluctu- 
ations are taken into account, the more the Neel state 
is stable, even in the isotropic limit. This is unexpected 
as the checkerboard antiferromagnet was expected to be- 
have like the kagome antiferromagnet. In the latter case, 
the reduction of the local magnetization to zero by quan- 
tum fluctuations and geometry, within a flnite range of 
couplings- has been shown using two distinct numerical 
methodslia. It was found that going from the triangular 
limit to the kagome limit, there is a singular point associ- 
ated to a critical couphng (r^ = J'c/ J ~ 0.9), where local 
magnetization vanishes. Here, we obtain the same con- 
clusion within the LSW approximation but not within 
the mean fleld treatment of the Dyson-Maleev approach. 
Wether this discrepency comes from the quantum fluctu- 
ation that are better incorporated in the latter formal- 
ism, or comes from the mean field treatment of the DM 
representation is not clear. jMevertheless, it is interesting 
to note that a recent workEj shows the same tendency, 
i.e a transition from an S" = ground state to a stable 
Neel state with increasing spin size S. If the present cal- 
culations treats correctly quantum fluctuations (which is 
not proven in this paper), the frontier pbetween the ex- 
pected VBS state and the Neel orderingt^l would then be 
at = 1. Nevertheless, it should be stressed that the 
present work does not test wether other classical ground 



states could be stable at the isotropic point and have the 
same or a different energy then the Neel ground state. 
Consequently, an intrinsic or accidental degeneracy of 
magnetic ground states for S" > 1 is not ruled out by 
this approach. 



V. THE S = 1/2 CASE 

As shown in the preceding sections, it seems that only 
for S = 1/2 can the ground state be non-magnetic. To 
shed light on the behavior of the S = 1/2 checkerboard 
antiferromagnet, it is possible to use a "fermionic" SU(n) 
generalization of the SU(2) group [Note that n > 2 does 
not correspond to a higher spin representation of SU(2). 
Therefore, a mean-field solution of the SU(n) hamilto- 
nian is expected to be much closer to t he SU (2) physics 
than the classical description of sec. 



II B . This for- 



malisiruJias been previously introduced by Marston and 
AfHeckEy and applied to the square antiferromagnet. The 
advantage of this descciption is that it works on all lat- 
tices, even frustratedEil or not bipartite, and that it is 
well suited for the checkerboard antiferromagnet as spin 
order never occurs in the exactly solvable large-n limited, 
as it is expected in the range Tc < r < 1 for n — 2. At 
large n, the hamiltonian can be written as 



H = 




H.c) 



(8) 



where Xij = l-'/'^jc-^Sa 
and where 
pressed by 

Si.Sj — (1/2) X]q,/3 '^iS/^'^Z/sSa 



— {J/^)cfa<^ja bond variables 

spin-spin interaction has been ex- 
use of electron operators through 
constant. The x 

variables may be seen as valence bond operatorsS as 
{n/ J)xtjXij is the number of valence bonds on the link 
{ij). Therefore, in this formalism, every ground state 
is described by a map of these complex scalar fields. 
Several authors have concfjTitratcd on the mean-field 
solutions = 0")l3E3o and in particular, it is 

possible pto apply here a theorem established by D. 
RokhsaiO. First it is clear that the checkerboard 
antiferromagnet is "dimerizeable with respect to J" as 
"it is possible to partition the network into disjoint 
pairs of sites such that (a) every site belongs to one 
and only one pair, and (b) for each pair (i,j) whe have 
Jij = J". In the range J > J' it is possible to apply 
his theorem which shows that piean-field ground states 
are (spin-Peierls or box)-phasesa (see Fig. ^l]). Keeping 
in mind that the case of interest is for n = 2, wee look 
for the proposed states that lower the energy of the 
real S — 1/2 case. It is straightforward to see that 
box-phases lower this energy with respect to J, as soon 
as these phases are regular tilings of the underlying 



square lattice by disconnected J-squares (see Fig. 11). 
Taking now into account the J' couplings, only one of 
these phases survives, where J' links are not frustrated. 
At this point, it means that if the SU(n) description of 




(a) 




FIG. 11: (a) General box phase on the J square lattice of the 
anisotrogic checkerboard antiferromagnet (as shown in Fig. 4 
of Ref . h4) . Bold lines without arrows correspond to non zero 
Xij while thin lines, plain or dashed, correspond to Xij = 0. 
Bold lines with arrows indicate that Xij = X^^^^* where x is 
real. The product of Xij around each box is negative, indi- 
cating a flux of n. These boxes are the large-n equivalent of a 
4-spins singlet on a square, (b) Picture of a general box phase 
for the 5 = 1/2 case. Bold squares correspond to 4-spins sin- 
glets and bold lines to 2-spins singlets. If J' interactions are 
within a box, they must be taken into account and are also in 
bold as they enter the calculation of the energy. The different 
energies for squares with crossings, squares without crossings 
and links are respectively E = —1.5, E — —2 and E — —0.75. 
(c) Picture of the valence bond solid for the 5 = 1/2 case, 
expected to be the ground state. This state is two times de- 
generate and involves only squares without crossings. 



the SU(2) case is relevant, it suggests that for S = 1/2, 
in the range < r < 1, the ground state should be two 
times degenerate (because there are two choices for the 
tiling), consisting in a valence bond solid of disconnected 
squares without crossings (in the range < r < rc it 
is clear that this formalism is not relevant as the local 
magnetisation is finite). The proposed ground state 
(see Fig. |l|(c)) could be tested by comparing the VBS 
energy per site E = —0.5 and the correlations fuBetions 
with results of exact diagonalisations for exampleEs. 



VI. DISCUSSION AND CONCLUSION 

To conclude. Linear Spin Wave analysis has shown that 
away from the isotropic checkerboard lattice limit, there 
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is a finite range of coupling for wich the Neel state is 
unstable, whatever the value of the spin. Strikingly, 
when using a self consistent decoupling within the Dyson- 
Maleev formalism, it is shown that the Neel state is 
never destabilized, excepted for the S = ^ case. The 
stability of the Neel state with respect to other magne- 
tizecLatates has not been tested in this work. Recent 
workEj supports that this state should be the ground 
state for large 5*. It is compatible with the present re- 
sults although it is at Vjariance with results obtained with 
the infinite-componentE3 antiferromagnetic model on the 
chekerboard lattice {d oo). A possible explanation for 
thisEZi is that for large S, thece should be quantum order 
by disorder (as supported iiJla) but this order by disorder 
comes with an energy scale of J/S rather than J. There- 
fore, in the classical limit, the temperature at which the 
correlation length grows exponentially goes to zero. In 
other words, the limits T — > and S' — > oo (or d — > oo) 
should not commute. For the special S — ^ case, the 
fermionic large-n limit suggests that the ground state is 
a valence bond solid of disconnected 4-spins singlets on 
squares without crossings. That conclusion is consistent 



with results of Refs. |T9| and |2^. 

The present work shows that the isotropic checker- 
board antiferromagnet should display behaviors quali- 
tatively different from the one of the kagome and the 
pyrochlore antiferromagnet. Moreover, quantum fluctu- 
ations have to be treated with care as they can drive the 
system from a purely quantum ground state to a Neel 
like ground state. In this frame, the structure of the 
singlet and triplet towers of state should be extremally 
5-dependent and exact diagonalization of finite clusters 
with different spin size would be of high interest, espe- 
cially to test wether the descrepency already appears be- 
tween the S = ^ and the S — 1 cases. 
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